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Abstract. We introduce and discuss discrete two-dimensional models 
for XY spin systems and screw dislocations in crystals. We prove that, 
as the lattice spacing e tends to zero, the relevant energies in these 
models behave like a free energy in the complex Ginzburg-Landau the- 
ory of superconductivity, justifying in a rigorous mathematical language 
the analogies between screw dislocations in crystals and vortices in su- 
perconductors. To this purpose, we introduce a notion of asymptotic 
variational equivalence between families of functionals in the framework 
of L-convergence. We then prove that, in several scaling regimes, the 
complex Ginzburg-Landau, the XY spin system and the screw disloca- 
tion energy functionals are variationally equivalent. Exploiting such an 
equivalence between dislocations and vortices, we can show new results 
concerning the asymptotic behavior of screw dislocations in the | loge| 2 
energetic regime. 

Keywords: Crystals, Discrete-to-continuum limits, Analysis of microstruc- 
ture, Topological singularities, Calculus of Variations. 

2000 Mathematics Subject Classification: 49J45, 74N05, 74N15, 74G70, 
74G65, 74C15, 74B15, 74B10. 



Contents 



1. Introduction 

2. Overview of the models 

2.1. Ginzburg Landau functionals 

2.2. The discrete lattice 

2.3. XY spin systems 

2.4. Screw dislocations 

3. The topological singularities and energy functionals 

3.1. The Jacobian 

3.2. Discrete dislocations 

3.3. Discrete vortices 

3.4. The energy functionals 

4. The variational equivalence argument 

4.1. The notion of variational equivalence 

4.2. Some consequences of the variational equivalence 

4.3. A first example of equivalent families 

5. Variational equivalence between QC S , Xy e and ST> £ 

5.1. Proof of xy e x gc e 

5.2. Proof of SV £ X Xy s 



m 
m 

7 
7 

m 

8 

m 

ED 
ED 
EH 
El 
EH 
E31 
El 
E2 
EH 
EH 



i 



2 



R. ALICANDRO, M. CICALESE, AND M. PONSIGLIONE 



5.3. Proof of QL E X SV E [23] 

6. New results for the asymptotic of SD e [21] 

6.1. The F-convergence result for QC £ in the | log e| 2 regime [25] 

6.2. New results for homogenizing dislocations in the | loge| 2 regime [26] 

7. Further extensions and conclusions [5H] 

7.1. The core radius approach [28] 

7.2. Three dimensional dislocations [29] 
Acknowledgments [3U] 
References [30] 



1. Introduction 

Since the pioneering papers of Berezinskii [5], Kosterlitz [25] and Kosterlitz and 
Thoulcss [26 , there has been a great effort in studying physical systems exploit- 
ing BKT-phasc transitions; i.e., phase transitions mediated by the formation of 
topological singularities of the order parameter. This type of phase transitions 
characterizes several physical phenomena such as superfluidity superconductivity 
and plasticity (see [23], [27], [21], [25], [31]), while vortices in superconductivity and 
XY spin systems, as well as screw dislocations in crystals, provide three paradig- 
matic examples of singularities. 

The phenomenological analogies shown by these apparently far physical systems 
have been pointed out many times in the physical community. In the language 
of statistical mechanics it has also been rigorously proven that these systems be- 
long to the same universality class all of them sharing a BKT-type phase transition. 
Roughly speaking it is known that above some temperature threshold, these systems 
undergo a phase transition to a disordered state in which the topological singulari- 
ties are unbound. Under that threshold the correlation length exponentially decays 
and the singularities bind together and interact through complex and mostly un- 
known phenomena involving many interacting scales. Such a complex behavior is 
the main reason why a detailed analysis of the ground states of these systems turns 
out to be a non trivial task. Moreover, the above description explains, to some ex- 
tent, why a qualitative approach to the study of the thermodynamic limit of these 
systems, such as the celebrated Ginzburg-Landau theory, has been so successfully 
exploited. 

In this paper we are concerned with the problem of describing some relevant 
properties of the ground states of these systems in the thermodynamic limit. We 
aim to provide a unifying mathematical point of view, based on a variational equiv- 
alence argument, to study the asymptotic behavior of the ground states of different 
models that share the same geometrical and topological qualitative features. More 
specifically our purpose is two- fold. On one hand, we want to reinterpret several 
known results about the asymptotic behavior of the ground states of such models, 
proving that the corresponding free energies are indeed equivalent from a variational 
point of view. On the other hand, taking advantage of this equivalence, we want 
to exploit some of the results currently proved only through a phenomenological 
Ginzburg-Landau analysis, to obtain new results in different contexts. 

As in [26], among the physical systems exploiting topological type phase tran- 
sitions we focus on two-dimensional systems and we choose two paradigmatic ex- 
amples: screw dislocations (SD) and XY spin systems. We will introduce two 
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basic discrete models, both constructed on eZ R f2, where f2 C M 2 is a bounded 
open set. Their order parameters are a unit vectorial spin field for the XY model: 
i £ el? n O i — > v(i) S M 2 such that \v(i)\ = 1, and a scalar displacement field for 
the SD model: i S eZ 2 n il i— > u(i) £ HL For a given configuration of spins or 
displacements, the energies of these systems are given by 

XYM-=\ E K*)-^)l 2 

i,jesZ 2 nQ: \i-j\=e 

and 

SD e (u):= l - J2 dist 2 (u(i)-u(j),Z). 

ij"SeZ 2 nf2: |i-j|=e 

With the energies written in this form, the coarse-graining analysis now amounts 
to study the limit, as e — > 0, of (some scaled version of) XY e and SD e . To this 
purpose, in the physical literature, it is customary to perform a so-called Ginzburg- 
Landau (GL) analysis (see [H] and [31] for an introduction to the subject and some 
applications). The main ansatz of this approach (based on heuristic scaling and 
symmetry type arguments) is to assume that some of the interesting features of the 
thermodynamic limit of the original functionals can be obtained by studying the 
limit, as e — > 0, of a family of so-called complex Ginzburg-Landau energies. These 
energies have as order parameter a vectorial field x £ £1 i— > w(x) £ K 2 and are 
defined as 

(1) GL E ( W ):= [ l|V™| 2 + i(l-M 2 ) 2 . 

This kind of functionals has been originally introduced as a phenomenological 
phase-field type free-energy of a superconductor, near the superconducting tran- 
sition, in absence of an external magnetic field. Here the order parameter w de- 
scribes how deep the system is into the superconducting phase and the scale e is 
proportional to the coherence length of the superconductor. 

The Ginzburg-Landau functionals have deserved a great attention by the physical 
community. The first rigorous mathematical approach to the limit as e — > of the 
solutions to the Euler-Lagrange equation of (JT|) has been made by Bethuel, Brezis 
and Helein in [7]. Since their paper a great effort has been done to study the 
asymptotic behaviour of minimizers of the Ginzburg-Landau energy both from the 
PDEs and the Calculus of Variations points of view. In particular in dimension 
two Jerrard and Soner in [25] (see also Alberti, Baldo and Orlandi [T] for the 
generalization to any dimension) have proved a L-convergence result for ^^p- 
In their analysis the relevant tool to track energy concentration is the asymptotic 
behavior of the Jacobians J(w £ ) of sequences (w e ) equi-bounded in energy. In 
particular they prove that, up to subsequences, J(w s ) converges to a finite sum of 
Dirac deltas whose support represents the vortex-like singularities of the limit field 
and that the L-limit is proportional to the number of such singularities. 

Only recently a similar analysis in the context of spin systems and dislocations 
has attracted much attention in the mathematical community and it has been 
carried on both in a continuous framework (see [T3], [2], [TS], [20]) and in a discrete 
setting (see 0, [3], [2], [3D]). As a further remark we underline that the asymptotic 
analysis of spin systems and discrete dislocations energy functionals is itself part 
of a wider interest in the discrete-to-continuum limits for more general models (see 
for example [5], [S] and [TT] (Chapter 11) for a review on this subject). 
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In [3] and in [30] a T-convergence result for Xy-spin systems and for the SD 
model is given in the | log e | scaling regime. Here the T-convergence analysis is 
performed with respect to the convergence of the Jacobians of a suitable affine 
interpolation of the spin variable v for the XY model, and with respect to the 
convergence of a suitable discrete notion of the curl of the strain field u for the SD 
model. Roughly speaking, gathering together the main results of these two papers, 
the following relations hold: 

r- hm = r- hm = r . lim . 

£^0 | log £ | £^0 |log£| e^O |loge| 

Motivated by this chain of equalities, we were led to ask whether one could pro- 
vide, in the framework of T-convergence, a unifying mathematical point of view to 
rigorously relate the asymptotic behavior of these models. Our purpose is to prove 
that the asymptotic equivalence of these models, in terms of T-convergence, can be 
push forward to any | log el' 1 scaling regime with h > 1. More precisely, we show 
that 

(2 ) r- hm = r- hm = r- lim ^ SD f\ 

e^o | log e\ h £^o|loge|' 1 e^o \loge\ h 

In this way we rigorously obtain the equivalence of some mean field models for 
vortices, spin systems and dislocations, according with experimental evidence (see 
|31j for a recent overview of the analogies between the mean fields in these mod- 
els). In particular, we obtain a rigorous justification to the GL analysis of the 
thermodynamic limits of the XY and SD models in the | log el' 1 energetic regimes. 

To prove ^ , we look for a relation between the order parameters of the different 
models, at proper mesoscopic scales, which allow us to compare the three families 
of energy functionals in the | logs I' 1 scaling regime for every h > 1. This has led 
us to introduce a notion of variational equivalence between families of functionals 
(see Definition [7]) . To explain the meaning of such a notion, let us suppose that we 
are given two families of energies (F £ ) and (G £ ) depending on a small parameter e 
standing for an interaction scale. Then we say that (F e ) and (G e ) are variationally 
equivalent if there holds that G e ^ F e and F e < G e . By G £ < F E we mean that, 
for any given family of order parameters (p E ) such that F E (p e ) < C, there exists 
another scale 8 S and a family of order parameters q e such that q e is closer and 
closer to p e and Gs E (q e ) < F £ (p e ) + 0(1). Roughly speaking, we are saying that 
the two variational models whose energies are given by F e and G e describe the 
same phenomena if looked at proper interaction scales and by suitably choosing 
the order parameters on these scales. The main feature of this notion, is that 
variational equivalent families of functionals share the same T-limit and share the 
so called equi-coercivity property (see Theorems [8] and 10 1. 



We remark that our notion departs from the concept of asymptotically equivalent 
functionals (at a certain order) introduced by Braides and Truskinovsky in |12j . 
Roughly speaking, according to their definition, given a > the functionals F e and 
G e are said to be equivalent at order a if % and % have the same T-limit. In 
particular our definition turns out to agree with the latter at order zero provided 
(F e ) and (G e ) are equi-coercive. On the other hand, to simplify matter, the purpose 
of the authors in [12] is to introduce a formalism to build up variational models 
that share the same T-limit up to a certain given scaling order, and then to study 
the properties that such a convergence enjoys with respect to a given family of 
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Figure 1. A screw dislocation "corresponds" to a vortex of spins. 



parameters specific of the considered theory. Our aim is instead to deduce T- 
convergence and compactness results for the family G e from the same results for 
the equivalent family F £ . 

In Theorem [15] we prove the equivalence betwe en t he families of functionals 

GL C \ ( XY e \ „„ i ,iTr 2 SD E 



(^), (j^pr) and (^f) (see also Theorem 



2:-! 



for h = 2). The way our 



variational equivalence is proved provides an interesting identification between the 
order parameters and the corresponding singularities of the different models. For 
instance, the identification underlying the equivalence between | and | 

relies on suitable interpolation procedures which allow us to pass from the discrete 
order parameter of the XY model to the continuous one of the GL model. Analo- 
gously the displacement field in the SD model is identified with the phase function 
of the XY order parameter (see also Remark [2| . These identifications of the order 
parameters clearly induce an identification of the corresponding singularities (see 
picture [T|). In the proof of Theorem 15 we have to make sure that the identifica- 
tion of the order parameters in the different models produces small perturbations 
in the corresponding free energy densities. This is easily checked far from the sin- 
gularities, while to control the error near the singularities we need to introduce 
suitable reparametrizations 5 e of the correlation length, i.e., we have to look at 
the models at suitable meso-scales 5 £ . Finally, during these identifications, we have 
also to control the distance, measured in a suitable topology, of the corresponding 
singularities. This analysis involves notions of geometric measure theory, and the 
arguments used in the proofs are close to those used to prove density of polyhedral 
boundaries in the space of integer currents in [T7] and also exploited in pQ . 

Taking advantage of this equivalence principle, we are able to export many of 
the known results in the theory of GL vortices to the framework of SD models. 
Indeed the T-convergence results in [23 in the | log £ | 2 energetic regime, together 
with ([2]) for h = 2, leads to new asymptotic results in the context of dislocations 
and spin systems when the number of defects grows logarithmically as e goes to 
zero (see Theorem 24 1 . 

The I log e\ 2 energetic regime has been already considered in the vectorial context 
of homogenizing edge dislocations in [TI5] , within a core radius approach, under the 
assumption that the dislocations have a minimal distance of the order of a suitable 
meso-scale. The T-convergence analysis done in |19j provides a macroscopic model 
for plasticity, in agreement with the phenomenological strain gradient theory for 
plasticity introduced in 1 18] . Moreover, the limit energy is compatible with the 
experimental evidence of the concentration of dislocations on lines, usually referred 
to as dislocation walls (we refer to [TS] for a variational model describing dislocation 
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patterns in crystals). In view of ([2]), we extend the T-convergence analysis done in 
|19j to our completely discrete setting, without any kinematic assumption on the 
mutual distance of the screw dislocations. In this way, we derive a strain gradient 
model for plasticity in the scalar setting of anti-planar elasticity, starting from a 
completely discrete and basic model of screw dislocations. 

Finally, let us mention that, as a byproduct of our equivalence result Q, taking 
into account the T-convergence result proved in [3D], we obtain (see Remark [18]) a 
new proof of compactness and T-convergence for two-dimensional GL functionals, 
independent of the proof given in J33] and in [35J . 

Our method can be clearly exploited in many interesting directions: first, one can 
investigate the equivalence of lower order terms for the models we have discussed 
so far, comparing the so called renormalized energy for vortices and dislocations 
within a T-convergence analysis, in the spirit of the theory of development by T- 
convergence introduced by Braides and Truvskinovsky in |12j . Moreover, one can 
consider the case of edge dislocations, or, more in general, the case of three dimen- 
sional models. Indeed, the results of this paper provide a first step in the effort of 
making a link between material dependent models for dislocations and phenomeno- 
logical Ginzburg-Landau approaches. We believe that our arguments could give 
efficient hints to build up material dependent Ginzburg-Landau energies, taking 
into account kinematic constraints and elasticity constants specific of the crystal. 
Moreover, exploiting our variational equivalence arguments in the three dimensional 
problem (e.g., in a cubic crystal) would bring new light on interesting mathematical 
questions regarding compactness properties and asymptotic behaviour of general- 
ized Ginzburg-Landau functionals, the target space being a three-dimensional torus, 
and the singularities being rectifiable currents with multiplicity in the group Z 3 (see 
Section 

The paper is organized as follows. In Section [2] we introduce the discrete mod- 
els for spin systems and screw dislocations, while in Section [3] we introduce the 
corresponding topological singularities. In Section [4] we describe our variational ar- 
gument, that will be used in Section [5]to prove the variational equivalence between 
GL £ , XY e and SD e models. Such an equivalence will be specialized in Section [6] 
in order to present new results in the asymptotic analysis of screw dislocations. In 
Section [7J we will comment the results achieved in this paper suggesting further 
extensions, concerning, for instance, the core radius approach to the singularities. 
Finally we will propose, in the case of three dimensional elasticity, a material de- 
pendent Ginzburg-Landau type model for dislocations in a cubic lattice. 

2. Overview of the models 

In this Section we briefly describe the models of Ginzburg-Landau vortices, of 
XY-spin systems and of screw dislocations. We will provide a detailed description 
of the latter in order to define the physical quantities involved in the model and 
needed to correctly describe the new results in the framework of screw dislocations 
contained in Section [6] 

For the time being f! C I 2 is a bounded open set with Lipschitz boundary, 
representing the domain of definition of the relevant fields in these models. For 
the sake of simplicity, we will also assume that f2 is star-shaped with respect to 
the origin. We stress that with some minor technical effort in our proofs, such 
assumption can be removed. 
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2.1. Ginzburg Landau functionals. Let us introduce the family of the so-called 
complex Ginzburg-Landau functionals GL e : VF 1,2 (f2; R 2 ) — > [0, +00), defined as 

(3) GL e {w):= ni\Vw\ 2 +±W(w) 

where W(x) :— (1 — \x\) 2 . Here w represents the order parameter of the model, 
describing how deep the material is in the superconductive phase, e is a length-scale 
parameter, usually referred to as the coherence length while GL e is the correspond- 
ing free energy of the system. 

Remark 1. We make this explicit choice for W = (1 — \x\) 2 because it simplifies 
some computation. The following standard hypotheses would suffice to perform 
our analysis: W G C(R 2 ) such that W{x) > 0, = S 1 and 

hmmf — > 0, hmmf , ]' > 0. 

|o!|-l (1 - \x\) 2 |xHoo \x\ 2 

2.2. The discrete lattice. Here we introduce the discrete objects and notations 
we will use in the sequel. 

For every positive e > 0, we set := eZ 2 nfi, representing the reference lattice. 
We will denote by fi* := £ 0° x 0° : \i — j\ = e, i < j} the class of nearest 

neighbors in Sl^ (where i < j means that i\ < ji for I G {1, 2}). The class of cells 
contained in is labelled by the set £l 2 :— {i G £1° : i + [0,e] 2 C £1}. Finally, we 
set 

n E -.= J |J J + [0,£] 2 
In the following, we will extend the use of these notations to any given open subset 

A of m 2 . 



2.3. XY spin systems. Here we recall the model of XY spin system following 
the approach in [3] (see also [M] for a general introduction to the model) . First we 
introduce the class of admissible fields, 

(4) AXy E := {v : S 1 }, 

where S 1 denotes the set of unit vectors in R 2 . The family of functionals XY e : 
AXy e -> R are defined by 

(5) XY e {v):=\ Hi)-v{j)\ 2 . 

(i.j')eni 

We refer the interested reader to [3] for the derivation of these energies by proper 
scaling of the XY energies written in the usual form 

where (a, b) denotes the scalar product between the vectors a and 6. 
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2.4. Screw dislocations. Here we introduce a basic discrete model for screw dis- 
locations, inspired by the approach introduced in [5] and revisited in 30J . The 
displacement, in this discrete anti-planar setting, is a function u : — > HL We 
denote the class of all admissible displacements by 

(6) ASV £ := {u : ft° -> M}. 

We focus here on linearized elasticity, and we consider the model case of nearest 
neighbors interactions, so that the discrete elastic energy corresponding to any 
displacement u, in absence of dislocations, is given by (we fix the shear modulus 

(7) Ef{u):=\ W*')-<i)| 2 . 

It is convenient to introduce also the notion of discrete gradient du, defined on the 
nearest neighbors (namely, the bonds of the lattice), by duij = u(j) — u(i), for 
every € 0*. With respect to the discrete gradient du, the energy ^ reads 

like 

(8) Ef{u):=\ Yl \ du ^- 

To introduce the dislocations in this framework, we adopt the point of view of 
the additive decomposition of the gradient of the displacement in an elastic part, 
the strain, and a plastic part, following the formalism of the discrete pre-existing 
strains as in [5] and |30j . More precisely, a pre-existing strain is a function (3 P 
representing the plastic part of the strain defined on pairs of nearest neighbors and 
valued in Z\b\. Here b represents the so called Burgers vector which is characteristic 
of the crystal. In principle |6| should be of order s, but up to a further re-scaling 
in the energy functionals, we can fix from now on |&J = 1. 

Here the idea is that the plastic strain (3 P does not store clastic energy and hence 
it has to be subtracted to the gradient of the displacement in order to obtain the 
so called elastic strain. In view of this additive decomposition Vw = /3 P + (3 e , we 
have that the elastic strain f3 e is not curl-free (in a suitable discrete sense). In 
Section [3] we will introduce the quantity /i := curl (3 e = —curl /3 P , that measures 
the degree of incompatibility of the elastic strain p e from being a discrete gradient, 
and represents the discrete screw dislocations in the crystal lattice. 

Summarizing, the elastic energy corresponding to the decomposition Vu = fP + 
(3 e is given by 

Ef(u):=\ £ \PU 2 = \ £ |d^--^.| 2 , 

and whenever the dislocation density curl [3 e = —curl (3 P is non zero, then the 
corresponding elastic energy is also non zero. 

Given a displacement u, we can minimize the corresponding elastic energy with 
respect to the plastic strain. Since by our kinematic assumption the plastic strain 
takes values in Z, it is clear that the optimal (3 P is obtained projecting du on Z. 
More precisely, let P : M. — > Z be the projection operator defined by 



(9) 



P(t) = argmin{|£ — s\,s € Z}, 
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with the convention that, if the argmin is not unique, then we choose that with 
minimal modulus. Then we have /3Jj = Pdu, in the sense that, for all G f2*, 
(Pu)i,j — P{dui,j)- In this way, we obtain the elastic energy functionals SD £ : 
ASV £ -> R defined by 

(10) SD e (u):=± J2 \(PZhi\ 2 = l E dtat a («(*)-«(i),Z). 

We notice that, defining u := eu as the physical displacement, we have 
e 2 5£> e («) = i £ dist 2 (u(z)-u(j),eZ), 

where the r.h.s. now reads as a term penalizing a misfit slip in the crystal, according 
with Pierls-Nabarro theories. 



Indeed, another way of understanding (10 1 is the following. In our anti-planar 
setting, the crystal can displace only in the vertical direction, and each vertical line 
of atoms has to displace rigidly. The discrete gradient du measures the difference 
of the displacement of two near lines of atoms. If du is of order s, the periodic 
structure of the crystal is unperturbed, and therefore the corresponding stored 
energy has to vanish. In other words, the nearest neighbors interactions have to be 
labeled in the deformed configuration, and not in the reference one. The rigorous 
way of formalizing this idea is given exactly by the projection procedure introduced 



by the operator P in ( 10 1 



Remark 2. Here we describe the heuristic argument to identify the XY and SD 
models just introduced. The correspondence between the displacement functions 
u £ G AST> £ and the XY fields v e G AXy e is given by identifying u £ with the phase 
function of v £ . More precisely, given a displacement u £ G AST> £ , the corresponding 
field v(u £ ) for the XY model is given by 

(11) v(u £ )(l) := e 27riu *W for every I G fi°. 

Viceversa, given v £ G AXy £ , the corresponding displacement is u(v £ ) :— ^-9 Ve) 
where 6 Vc G [0, 2ii) is defined by the identity v £ — e . Note that the (arbitrary) 
choice of a precise representative of the phase 9 Ve of v £ does not affect the elastic 
energy corresponding to u(v £ ). Indeed, we have 

27rdist(u £ (i) - u £ (j),Z) = d s i(v(u £ )(i) - v(u £ )(j)), 

where d s i denotes the geodesic distance on S 1 . 

Finally, we notice that by Taylor expansion we have 2ir{u s (i)—u e {j)) w v(u £ )(i) — 
v(u £ )(j), whenever u £ (i) — u £ (j) is small. Therefore, we expect the identification 
between the fields to produces small perturbations in the corresponding energy den- 
sities (up to a pre-factor Att 2 ) far from the singularities, and this will be formalized 
in the sequel. 



3. The topological singularities and energy functionals 



As discovered by Jerrard in [5T] (see also [H]), the Jacobian of the order pa- 
rameter is the relevant geometric object carrying energetic informations for the 
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Ginzburg-Landau functionals. In the same spirit, in this paper we introduce suit- 
able discrete notions for the topological singularities of the models we have in- 
troduced, and we consider them as the meaningful variables of the corresponding 
discrete energy functionals. 

For the SD model, the topological singularity is given by the discrete dislocation, 
defined through a discrete version of the curl of the elastic strain, while for XY 
spin systems, the natural notion of topological singularity is that of discrete vortex, 
that we will introduce in the sequel. 

We will follow the formalism introduced in [5] (see also [3D]) representing the 
singularities by a discrete function a whose values on the cells of the lattice represent 
the topological degree of the singularity. Moreover, we will identify this function a 
with a measure concentrated on the center of the cells of the lattice. 

As we will see in the sequel, we can always pass from a discrete representation 
of the singularities to a continuous one (and viceversa) by mean of interpolations 
procedures (and projection on finite elements space respectively). This will also 
be consistent with the topology we are going to use to perform the T-convergence 
analysis: the distance, measured in such a topology, between the continuous rep- 
resentation of the singularities and its discrete counterpart will turn out to be 
vanishing for sequences of order parameters with bounded energy (see Proposition 



161 



3.1. The Jacobian. Given w € i/ 1 (fi;]R 2 ), the Jacobian of w is the L 1 function 
defined by 

Jw :— detVv. 

Let us denote by Cy(Q) the space of Lipschitz continuous functions on f2 with 
compact support, and by X its dual. The dual norm of X will be denoted by || • ||. 
For every w £ H 1 ^; M 2 ), we can consider Jw as an element of X by setting 

<Jw,ip>:= / Jwtpdx for every ip € C° ,:L (f2). 
Jn 

Note that Jw can be written (in the sense of distributions) in a divergence form as 

(12) Jw = div (wi(w 2 ) X2 , -w 1 (w 2 ) Xl ), 

or cquivalcntly, in the form Jw = curl (wi Vu>2). By a density argument, we deduce 
that for every ip € C° ,:L (fi), 

(13) <Jw,tp>=- w 1 (w 2 ) X2 (p Xl - w 1 (w 2 ) Xl (p X2 dx. 



Note that the right hand side of ([13]) is well defined also when 

we W 1A (n;R 2 )nL°°(n;R 2 ), 



and for such a function, we will take ( 13 1 as the definition of Jw as an element of 
X. 

Finally, for later use we notice that for every v := (vi, v 2 ), w := (w%, w 2 ) belong- 
ing to H 1 ^^ 2 ) (or, as well, to W 1 ' 1 (Q;M 2 ) n L°°(0;IR 2 )), we have 

(14) Jv - Jw = ^{J( v i - »i, v 2 + w 2 ) - J(v 2 - w 2 , vi + wi)). 



By (13 1 and (14 1 we immediately deduce the following lemma. 
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Lemma 3. Let v n and w n be two sequences in 7? 1 (O; IR. 2 ) such that 

IK - w n || 2 (||Vu„||2 + ||Vu?„|| 2 ) -> 0. 
Then \\Jv n — Jw n \\ — > in X. 

3.2. Discrete dislocations. Following the formalism introduced in [5] (see also 
|30j). given a function £ : Q,\ — > K (playing the role of a pre-existing strain), we 
introduce its discrete curl d£ : £l 2 — * R, defined for every i 6 fi 2 . by 

(15) d£(i) := 6i,i+(0,e) + Ci+(Q,e),i+(e,e) - £i+(e,0),i+(e,e) ~ £i,i+(e,0)- 

Given an admissible displacement u : 0° — > R, we recall that we can decompose 
du in its elastic and plastic part (that is optimal in energy), by setting (3? = Pdu, 
/3£ = du — ft?, where P is defined in Q. 

We are now in a position to introduce the discrete dislocation function a u : f2 2 — ► 
{-1, 0, 1}, defined by 

a u (i) := d/3®(i) for every i £ fi 2 . 

It will be convenient also to represent a u as a sum of Dirac masses with weights in 
{— 1, 0, 1} and supported on the centers of the squares where a ^ 0, setting 

(16) Hu := a«(*)^+i(e, e )- 

Remark 4. Notice that, in view of the very definition of P, we have that (3^ £ 
(—1/2, 1/2], therefore by definition (15l the discrete dislocation a u takes values in 
{ — 1, 0, 1}. We deduce that, in our model, only singular dislocations can be present 
in a cell. 

3.3. Discrete vortices. Given an admissible field v £ AXy e , the associated dis- 
crete vorticity j v : f2 2 — > Z is defined by 

where 9{v) £ [0, 2tt) is the phase of v defined by the relation v = e 10 ^. Moreover, 
we introduce the measure fi v defined by 

( 18 ) M« : = 7o(*)^»+i( e , e )- 

Notice that, as for the screw dislocations, in a cell we can have only singular vortices. 

Remark 5. Given v £ AXy e , we can introduce a function v : il e — > R 2 that 
coincides with v on fl® n £l e and such that Jv = k(A v . Indeed, consider the function 
6(v) defined on each segment with £ ill by 

6(v)(i + s(j - »)) = 9(v)(i) + 2irs (^ {v) ) id for every s £ [0, 1], 

where u(v) := ^-0{v). Now, extend in each cell i+ [0, e] 2 (with i £ f2 2 ), making 
it zero-homogeneous with respect to the center of the cell. Finally, on each cell we 
set v := e l ^ v \ 

A straightforward computation leads to the equality Jv = ~k\x v . This argument 
shows that the vorticity function a v represents a discrete version of the Jacobian, 
and seems a very natural object in this context. 
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Finally, for latter use, we notice that it can be easily proved (for instance by the 
identity Jv — curl (5i VS2)) that for each e-square Qi with i G f2 2 we have 

If d 

(19) iJL v {Qi) = -l v 1 —v 2 ds. 

n JdQ, " s 

Remark 6. We notice here that a very easy estimate leads to the following bound 
for the total variation of the topological singularities 

\Jw s \{Q)<CGL s (w e ), \ti Ve \(Q)<CXY e (v e ), |/i u J(fi) < C SD £ (u e ). 

3.4. The energy functionals. In this paragraph we will formally introduce the 
(rescaled) energy functionals corresponding to the XY spin system, the screw dis- 
locations and the Ginzburg-Landau models written in terms of the corresponding 
singularities. 

In order to rewrite the energy functionals in these new variables, we minimize the 
free energies among all field quantities compatible with the prescribed singularity. 
For instance, in the SD model, we will fix fj, := curl /3 e , and minimize the clastic 
energy among all u compatible with /1, i.e., with \x u — [i. The corresponding energy 
functional can be thought of as the energy stored in the crystal, for a certain given 
dislocation density. 

We will consider the energetic regimes of order | log el' 1 , where h is any fixed 
positive real number. Following the convention according to which the infimum of 
the empty set is +00, we define the Ginzburg-Landau energy functional QC E : X — > 
[0,+oc], where X = ((^(fi))*, as 

(20) QLM := inf {gL £ (w),w G H\n^) : - { = M } . 

Let us pass to the energy functionals corresponding to XY spin systems. Using 

the notations introduced in the Section [3| the energy functionals Xy e : X — > 
[0, +00] are defined by 

(21) Xy e W := inf {xY e (v), v e AXy e : ^^^j^ = /i} . 

Finally the energy functionals corresponding to the screw dislocations model 
ST> E : X — > [0, +00] are defined by 

( 22 ) SV ^) := inf {SD^UI u G ASV e : = M } , 

where the pref actor 47r 2 is just a normalization factor which guarantees that the 
family (SV e ) asymptotically behaves as Xy e and QC e . 

4. The variational equivalence argument 

In this section we will introduce a notion of equivalence between families of 
functionals defined on a metric space {X, d) , depending on a small parameter e. 
Such a notion turns out to be efficient to compare different variational models 
which share the same asymptotic behavior as e goes to zero. 
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4.1. The notion of variational equivalence. Let (F e ) and (G e ) be two families 
of functionals from X to R U {00} depending on the parameter e G M. + U {0}. 

Definition 7. We set (G £ ) < (F e ) if there exists a continuous increasing function 
e 1— > S £ , with S = 0, such that the following holds. 

For every e„ — * 0, and (p n ) C X such that F £n (p n ) < C, there exists a family 
(q n ) C X such that 

i) limsup„(G 5en (g n ) - F e „(p„)) < 0; 

ii) Either (p„) and (g„) are unbounded or d{p n , q n ) — > as ri — * +00. 

We set (F e ) ~ (G e ), and we say that (F £ ) and (G e ) are variationally equivalent 
(for e - 0) if (F e ) ^ (G e ) and (G E ) r< (^)- 

Note that the relation ^ just introduced is transitive, i.e., if i 7 ^ 1 ^ i 7 ^ and 
F e 2 ^ F E 3 then F} ^ F e 3 . Moreover the relation ~ is an equivalence relation between 
families of functionals {(F e ),F e : X — + M.}, i.e., the following properties hold 

Reflexivity: (F E ) ~ (F E ); 

Symmetry: (F £ ) ~ (G £ ) implies (G e ) ~ (F e ); 

Transitivity: (F s ) ~ (G e ) and (G e ) ~ (ff e ) imply (G e ) ~ (ff e ). 



4.2. Some consequences of the variational equivalence. Let us consider now 
a first important relation between the notion of variational equivalence and that 
of L-convergence (for the definition and the main properties of L-convergence we 
refer the reader to [TO] and [I§] ) . To this purpose we recall the following definition 
of equi-coercivity: a family of functionals (F £ ) is said to be equi-coercive if, given 
£„ — > and (p n ) such that F Sn (p n ) < C, then (p n ) is relatively compact in X. 

Theorem 8. Let (F e ) be a family of equi-coercive functionals I '-converging to some 
functional H in X . Then (F e ) is variationally equivalent to the constant sequence 
H £ = H. 

Proof. We begin proving F e ■< H by a contradiction argument. Assume by contra- 
diction that the relation F e < H does not hold; then there exist sequences e n — > 0, 
(p n ) C X with H(p n ) < C, and r > such that 

(23) inf F e Jx)>H( Pn )+r. 

xeB r ( Pn ) 

On the other hand, by the equi-coercivity of F £ we easily deduce the coercivity 
of H , so that we can assume without loss of generality that p n converges to some 
By the L-limsup inequality there exists a sequence (z n ) C X such that 

(24) z n ^p, F En (z n ) ^ H(p)< limmfH{ Pn ). 



Equation (23 1 together with (24 1 provides a contradiction. 

Let us pass to the proof of the opposite inequality H -< F s . Assume by contra- 
diction that there exist sequences e n — > 0, (p„) C X with F £n (p n ) < C and r > 
such that 



(25) 



inf H(x) > F e (p n ) + r. 

xGB r (p n ) 
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By the equi-coercivity property of F £ we may assume that p n — > p in X, and by 
the L-liminf inequality we have that, for n big enough 

inf H(x) < H(p) < liminf F £ (p n ), 

xGB r (p„) 



which, together with (25) provides a contradiction. □ 

Remark 9. Note that in Theorem [8] the equicoercivity assumption plays a funda- 
mental role. Indeed consider the space I2 of sequences {a^} with ^ i af finite, and 
consider the functionals F £ : R — > K defined by 

F E ({ ai }):= J2 a *> 

where [1/e] denotes the integer part of 1/e. Clearly F e L-converge to the functional 
H({cii}) :— ||{a.,}|||, but F £ is not variationally equivalent to H . 

The next theorem, together with Theorem [8] clarifies the relation between the 
notion of L-convergence and the notion of variational equivalence. 

Theorem 10. Let {F £ ) and (G £ ) be variationally equivalent. Then the following 
properties hold. 

1) F £ are equi-coercive if and only if G £ are equi- coercive; 

2) F e L '-converge to some functional H in X if and only if G £ Y -converge to 
H. 

Proof. To prove property 1) assume for instance that G £ are equi-coercive and let 
us prove that so are also F £ . To this aim let e 1— > S £ be the map given in the 
definition of G £ ■< F £ . Let e„ — > and (p n ) C X be such that F £n (p n ) < C. By 
Definition [7] there exists a sequence (q n ) C X such that limsup„ Gs En (q n ) < C. By 
the equi-coercivity property of G £ we deduce that (up to a subsequence) q n — > x 
for some x G X. By property ii) of Definition |T] we deduce that d(p n , q n ) — > and 
therefore that p„-ti. 

Let us pass to property 2). Assume for instance that F £ T-converge to H and 
let us prove that also G £ T-converge to H. To this purpose, let e n — > and let p n 
be such that e n = 5„ n . In order to prove the L-limsup inequality let p£l, and let 
(p n ) be a recovery sequence for F Pn ; i.e., p n — > p and F Pri (p n ) — > F(p). Consider 
the sequences (g n ) C A" given by Definition [7j By property ii) of Definition [7] we 
have that q n — > p, and by property i) we have that 

limsupG e?i (<j„) < lira F Pn (p n ) = H(p). 

n 

In order to prove the L-liminf inequality, let now e —> S £ be the map given in the 
definition of F £ ^< G £ and let p n — > p. By Definition [7] there exists a sequence 
Qn ~ * P such that 

iJ(p) < liminf F s (q n ) < liminf G £n (p n ), 

n n 

that proves the L-liminf inequality. □ 

Remark 11. Note that if F £ and G £ are variationally equivalent, by Theorem [T0| 
we deduce that they share the same L-limits even when they are not equi-coercive. 

We also notice that in the class of equi-coercive functionals admitting a L-limit, 
our definition of equivalence coincides with that given (at order zero) in [12 . 
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Remark 12. Even if Dcfinition[7]is given in a metric context, it can be generalized 
to vectorial topological spaces. We chose the metric framework because it is usu- 
ally general enough for practical applications, the metric d being either the metric 
inducing the weak topology of a ball in a Banach space, or (in view of a compact 
embedding) the distance induced by the norm of a Banach space. 

4.3. A first example of equivalent families. For all positive s > consider the 
Ginzburg-Landau functionals GL S E : W 1,2 (£l; K 2 ) — > [0,+oo), defined as 



(26) 



GL s e (w) 



;W(W). 



Note that for s = 1 the functionals GL s e coincides with the functional GL £ defined 
in Consider also the functional QC s e : X — ► [0, +oo] defined as in (20) with GL e 
replaced by GL%. 

Proposition 13. For every < s\ < S2, the families of functionals QC^ 1 and QC S £ 2 
are variationally equivalent, according to Definition^ 

Proof. Let si < s 2 . Since QC^ 1 < G£° 2 , we immediately deduce QC.% 1 < GC S £ 2 . In 
order to prove the opposite relation, consider the following change of variables 



Following Definition [7j let e n 
gCl 1 (p n ) < C. Moreover, set 



e^5 e := e 
— > as n 

|i°g^„l fe - 1 



-oo and let (p n ) C X be such that 



._ Pn 

|log£ rl |' 1 - 1 ' ? " i ■ 

By the fact that t n — > 1 as n — > +oo it immediately follows that either p n and q n 
are both unbounded or d(p n , q n ) 
Let (w n ) C H 1 ^ 2 ; 

Then, since by construction j |h-i Jw n = q n , we have 



0. 

I 2 ) be such that 



-oo. 



limsup(0£f (q n )-GC s e l(Pn)) 



< limsup 

n 

lim sup ( 

|log(e 



logoV 



|loge n |' i 



I log e„| 



si 



\Vw n \ 2 + ^W{w n ) = 



□ 



In the following we will need a slight variant of Proposition 13 where the pre- 
factor s in front of the potential term W is a suitable function of e. 



Proposition 14. Let e i— ► s e oe an increasing functions from 



to 



such that 



I logS e | 

|loge| 



as e — > 0. Then the functionals QL e are variationally equivalent to 



QC S S " , according to Definition^ 
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Proof. The proof follows the lines of the proof of Proposition [13] setting now S e :— 



5. Variational equivalence between QC £l Xy £ and SV £ 

In this Section we prove the main result of the paper: the energies corresponding 
to XY spin systems, to screw dislocations and to the Ginzburg-Landau model are 
variationally equivalent in the sense of Definition [7] We prove this result for any 
scaling of the energies of order | log e\ h with h > 1, the most relevant cases being 
h = 1 and h = 2. 

In order to define all the energy functionals in the same space, we consider the 
Banach space (X, || • ||) defined as the dual of Lipschitz functions with compact 
support in f2, endowed with the dual norm. This space is the natural space for the 
topological singularities, that are the relevant quantities in all the investigated mod- 
els. The functionals GC £ , Xy e and SV e from X to [0,+oo] are defined rigorously 



in (201, (21) and (22 1 respectively. 



Theorem 15. The functionals Xy e , ST> e and QC e are variationally equivalent 
according to Definition^ 

We will prove in three steps that Xy e -< QC e -< SV e ■< Xy e . This relations 
imply the equivalence of all the functionals by the transitivity property of the 



order relation F -< G. Before the proof of Theorem 15 we state and prove a 



proposition that will be used to identify the discrete description of the singularities 



{i.e., the discrete measures fj, u and /i v defined in ( 16 1 and ( 18 1 , respectively) with the 
diffuse Jacobians of some suitable interpolation of the corresponding fields. To this 
purpose, for every positive e > 0, consider the triangulation T e :— {T ± (i), i G el?} 
whose triangles arc of the type 

T~(i) := {(xi,x 2 ) ■ h < xi < ii + e, i 2 < x 2 < i2 + xi}; 
(27) T+(i) := (ii,ii+e) x (i 2l i 2 + e) \ Tr(i). 



To any v e in AXy e , we can associate the continuous field w(v E ) in iJ 1 (f2 e ; R 2 ) given 
by the piece- wise linear interpolation of v £ on the triangles of T £ contained in Cl e . 
Note that the Jacobian of w(v £ ) is piecewise constant and that it can be extended 
by zero to a L 1 function in fi. Therefore J(w(v e )) can be seen as an element of X. 

Since we will need to localize the energy functionals to subsets of fi, we introduce 
the following notation. Given A C O, we will denote by GL e (-,A) the restriction 
of the GL E energy density to A, and we denote by XY E (-,A) and SD e (-,A) the 
restriction of the corresponding energy densities to the nearest neighbors contained 
in A. Finally, for every given positive <5 we set 

(28) Isity '■= {x € Q '■ min min \xi — yA > 5}. 

yedn ie{i,2} 

Proposition 16. Let v £ £ AXy e be a sequence with XY £ (v £ ) < C\loge\ h . Then 
||^- M J|^0 ase^O. 
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Proof. Let S e := (e| loge| denote by S e the projection of S e on eZ, and set 



m :={(x,y) eR 2 :ye6 E Z}, 



= {(x,y) £ M 2 : x £ 5 £ Z}, 

= (R}uR})n i s - e (Q). 



If we denote by Rg (s) the class of nearest neighbors in Rg + (s, s), by the Mean 
Value Theorem it is easy to prove that there exists sg £ e{0, 1, . . . , ^} such that 

XY £ {v £ )> E \\ve(l)-v e {k)\ 2 



(29) 



e ^ 2' 



By ( 29 1 we deduce that 



max |We(i) — v e(k) I — > 



as e — > 0. 



Let us denote by Q l s ~ the cubes contained in Ig (CI) + (sg , sg ) of the type 

Qg e ■= S e (i + [0, l] 2 ) + (sg e ,sg e ), % £ Z 2 , 

and by Ag their union. 

Let <p £ C®' 1 ^) with norm less then one, and denote by <p the locally constant 
function that on each (5 e -square Q l - coincides with ip on the center of Q % - . Then 
we have 



(30) |<I^M_ M> | 



< 



\Jw(v £ )\dx+\(j, Ve \(Cl e \Ag) max \ip\ 
Qe\A 5 - I n s \A s - e 



E 

01- CA t 



\J(w(v e ))\dx + \n Ve \(Q} e )\ \<p-<p 
Jw(v e ) 



E 



dx - fi Ve (Qg ) (p. 



The first two addends of the right hand side of ( 30 1 are vanishing (uniformly with 
respect to <p belonging to C®' (CI) and with norm less then one), since maxj^^ \ip\ 

and \<p — <p\ are bounded by CS e , and (see Remark ^ 

\Jw(v e )\dx+\^ Vs \(Q e ) < C|loge|\ 



Therefore it remains to prove that also the third addend in ( 30 1 is vanishing 



uniformly with respect to ip. To this purpose, let v e be defined as in Remark 
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|5j on the boundaries of all Q l g . Since J(w(v £ )) = curl (w(v £ ))iW(w(v £ ))2 and 
Jv e = curl (ti e )iV(u e )2 (see also of (19l), we have 

(31) / Jw ^-dx-^ Ve (Q i 5 -) 



1 

n JdQ 
1 



d d 
(w(v £ )) 1 — (w(v £ )) 2 - (v e )i — {v e ) 2 



d 



2 _ (w(v e )-v e ) 1 — (w(v e ) + v e ) 2 



(w(v e )-v £ ) 2 Q s (w(v e )+v e ) 1 . 



By (29 1 and (31 1 we conclude 



e(7. ^-^uW^e/ . 



k(w e ) - v £ \\w(v £ y\ 



< c 



W H«e)-£ £ | 2 W \w(v e )f 

<c Yl Mi)-v e {k)\ a <c^ 



logs | 



where in the last but one inequality we have used that \w(v £ ) — v £ \ is controlled 
on the segment [l,k] by \v £ (l) — v £ (k)\. 

Remark 17. Note that, given w(v £ ), we can always extend it to a function w(v £ ) 
defined in the whole such that 

\\Vw(v £ )\\ L 2 (n . R 2x2 } < C||Vw(w e )|| L 2 (ae . R 2x2 ) . 



As a consequence, if w(v £ ) is with finite energy as in Proposition 16 then it is 
easy to see that || Jw(v £ ) — Jw(v £ )\\ — > as e — » 0. Therefore, the conclusion of 
Proposition 16 still holds with Jw(v £ ) in place of Jw(v £ ). 



We are now in a position to prove Theorem |15| By simplicity of notation in 
what follows we will replace e n by e. 

5.1. Proof of Xy e < QC £ . In order to follow Dcfinition[7]we have first to define the 
function e h- > S £ . A convenient choice for this purpose is to set 5 £ := (e\ loge\ h+1 ), 
and S £ := X £ 5 £l where a suitable factor A E — > 1 will be chosen in the following. Let 
(fi £ ) C X be a sequence such that G£ £ (fi £ ) < C. By the very definition of QC £ 
there exists (w £ ) C i? 1 (fi;IR 2 ) such that 



Jw F 



I log el' 1 1 7T 



= fi £ and 



|loge| 



: GL £ {w £ ) - QC £ (n £ ) -> as e -► 0. 



let / 5 ~ (f2) be defined as in ( p8| (with <5 replaced by <J e ), and consider the nets Bh , 
and Rg^ defined by 
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R'l :={{x,y) eR 2 : y e S £ Z}, 
R} := {(x,y) S R 2 : x e S £ Z}, 
R s ~ e := (4u]J|)nJ 4 ((]). 



By the Mean Value Theorem, it is easy to prove that for every 8 £ we can find 
Sg e (0, S £ ) 2 such that 



1. ,,9 1 , s f 1, ,9 1 

..' \£ I TJ/Y„.. \ I Y7„.. K 



(32) 5 e y -K| ^ + ^WK) < J -\Vw e \ < + ^W(w e ), 

where v/ E (x) denotes the tangential derivative of w £ \—(Rs e + sg ) in x (that is well 
defined by standard slicing arguments). 

Claim: let I be a segment with length larger than e and let ip £ iJ 1 (L;M 2 ), 
then 

(33) J^ + ^ W(V) >C ^MZ^ . 

We now use the Claim, that we will prove later, in order to conclude the proof of 

xy £ ■< QC £ . 

Set v £ and v e from Rg to R 2 as follows 



[l):=w £ {l + s s - c ), v e (l) -r|4|r for every I e 5 £ I? n J & (f>). 



Note that by the Claim ( 33 ) and by the choice of S £ we immediately deduce that 



\v £ \ — 1 — > uniformly as e — > 0. 

Here we use the assumption that is star-shaped with respect to the origin. Let d 
be the distance between and dil. We set 



(34) 



d 



d-2S E ' 



so that C \ £ I 2 g (O) and 0° C A e (<5 e Z 2 n /j (f2)) . We are in a position to 
introduce the function 

(35) v £ (l) :=%(—) for every I 6 . 
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By Jensen inequality and in view of ( 32 ) we have 
(36) XY Se (v e ) = \ J2 \v E {i)-v e (j)\ 2 



(t, 3 )enl 



(i+o(i)) Yl m*)-Uj)\ 2 



;\we(i + Sg ) -w e (j + s s ~ )\* 



< 4(l + o(l)) 



1 



i |2 



(Rs E +ss.) 
1 



< (l + o(l)) / -\Vw E \ 2 + ^W(w e ) = (l + o(l))GL s (w £ ), 



where we recall that (J 5 - (fi))*- denotes the pairs (i, j) € S £ Z 2 n (fi) with i < j 
and |i — j| = 1. 

We are in a position to introduce the sequence of variables (r) e ) C X for the 
Xy& c functional, satisfying properties i) and ii) of Definition [7j 



llog^l"- 1 ^" 

By (36), since p^^j — * 1, we deduce that property i) of Definition j^FJ is satisfied. 

Now we will prove that ||/x e — t £ rj £ \\ — > for some £ e — > 1, that will ensure 
property ii) of Definition [7] To this purpose, set 



te 



llogel"- 1 ' so that /i E - i £ ?y e — | loger |h-i 

In view of Proposition [16] to conclude it is enough to show that 
(37) \\ 0& \\h-i pw e - Jw(v e )\\ ->0 ase^O. 



Note that (see Remark 17 1, we can always extend w(v E ) to (and we will still 
denote this extension by w(v £ )) such that ||Vto(w s )||2 < C|loge|' i . Therefore, by 

(since we also have || Vw e ||| < C|loge|' 1 ), in order to 
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Lemm a |3| and Remark 
prove ( |37[ ) it is enough to check that 

| loge| h ||io(t; e ) - w e \\2 — * as e — * 0. 

Let J2,s e (f2) be defined as in ( p8| with <5 e replaced by 25 e (so that all the functions 
we have just introduced are defined on l2S e (Q))- Using that |f2 \ l26 c \ < CS £ and 
that the potential W in the GL £ functionals controls the L 2 norm of w £ , it is easy 
to prove that 

lloge^IKO - w e ||| 2( n\/ 2Sc (n);a=) °- 
Therefore, we will estimate the L 2 norm only on I 2 s e (ty- 
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By triangular inequality we have 
(38) f Hv £ )-w £ \ 2 <C'f \w(v £ )-w(v £ )\ 2 + \w(v £ )~w(v £ )\ 2 



+ C / \w(v £ ) - w £ (x + sg )r + \w £ (x + 8 B -)- w £ \ z . 
We can easily estimate the first and the last addend in the right-hand side of 



dx 



(|38| as follows 

(\w(v £ ) - w(v £ )\ 2 + \w £ (x + s s - e ) ~ w £ \ 2 ^ 

<C8 2 [ \\7w £ \ 2 dx < CS £ 2 \loge 
Jn 



h 

vw E \ ax vw) E |iugfc| 

n 

Let us pass to the second addend. For any £ (Ig (O)) j~ , set L i: j := sg 



Moreover set rriij :— ma,x t< zL i j \\w £ (t)\ — 1 . Then by ( 32 1 and by the Claim 



(33 1 we get 



f \w{v £ )-w{i £ )\ 2 <C% 2 



<CS £ 2 s Yl [ \\<\ 2 + ^W(w £ )<C5 £ e\loge\ h . 



Let us pass to estimate the third addend in the right hand side of (38 1. Let 
i £ If (fi) := {i € 6 £ Z 2 : i + [0, S £ ] 2 C Ig (Q)} and let f t : (0, 5 £ ) 2 -> M 2 be defined 
by 

fi(x) := w(v £ (x + i)) - w £ (x + + i). 
Therefore fi(0) — 0, and for every x — {x\, X2) G (0, S £ ) 2 we have 

\m\ 2 <cs £ (^J* E \^j l ( tl x 2 )\ 2 dt+ \^j t (o,t)\ 2 dt^ . 

Integrating with respect to x in (0, S £ ) 2 we have 

\h{x)\ 2 dx 

(OA) 2 



and hence 



<C8 £ [5 £ \Vfi(x)\Ux + 6 e / \—M0,t)\'dt 

'(0/ e ) 2 JO dx 2 , 



[ \w(v £ ) - w £ (- + s s J\ 2 dx < V / \Mx)\ 2 dx< 



(v £ )\ 2 



C(C / \Vw(v £ )\ 2 + \Vw £ (- + s g )\ 2 + S £ / \w'(i 

+ K(- + ^)| 2 ) <4 2 |io g£ | f 
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We conclude by proving the Claim (33 1. Let t € L be such that 
m := max||v?| - 1| = \\ip(t)\ - 1|, 

and let (a, b) C L be the maximal interval containing t such that ||y(s)| — 1| > 
for all s € (a, 6). Then we have 

(39) J l \W\ 2 + ^W{<p) > J l ^W(<p) > C l - 

if \b—a\ > e/2 we are done; otherwise either we have ||y>(o)| — 1| = ^ or — 1| = 

~. Then it is very easy to see that 

(40) j L \^' + ^WW)>j L \\9f>C^>C^, 

and this concludes the proof of the claim. □ 



5.2. Proof of SV E <Xy e . As in the proof of Xy E ^ ££ £ , we set (5 £ := e\ loge| h+1 , 
S e = \ E 8 E with \ E := - d ^ , where d is the distance between and dfl. Moreover, 

we denote by S E = e[^f], where [■] denotes the integer part. 

Let(/x £ ) C X be such that Xy e (ji e ) < C. By the very definition of Xy e there 
exists (v E ) C AXy e such that 

| log£ |ft-i ^ = ^, j^^xY E {v E ) - xy e {n e ) ^ as £ » 0. 

Let Ig (f2) be defined as in ( p8| (with (5 replaced by S E ), and set 

(41) fl* :={(ei,S e j),i,j eZ}, 

(42) i^:={(4* J ej),*,ieZ} > 

(43) % e :=(4 c Ui?yn/ 4 (0). 

Therefore, there exists sg S e{0, 1, . . . , ^} such that, denoting by R 1 ^ (s) the class 
of nearest neighbors (I, k) € + (s, s), we have 

(44) XY E (v E )> E hve(l)-v e (k)\ 2 



>- E Jk(0"-e(fc)| 2 . 

£ 2 

Let E {v E ) be a determination of the phase of v e , defined by the identity v E (l) — 
e ie e (v e ){l) f or ever y i e q0_ By (j44j) 5 smce XY E (v E ) < C\loge\ h , we immediately 
deduce that 

sup{|« e (0-i7 e (fc)|, (l,k) £ R} e {s s - s )}^0 ase^O, 
so that by Taylor expansion we have 
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2ir dist 



Q-e £ (v £ ){i) ~ ±-e e (v £ )(k) 

= e e (v e )(l) - 6 s (v e )(k) = (1 + o(l))|» 6 (0 - v £ (k)\. 



Let us set Ag :— I 2 g (f2) + (sg ,sg), and let (Ag ). be the class of <5 e -nearest 
neighbors in Ag . By Jensen inequality, in view also of ( 44 1 , we deduce 



(45) 



4^ 2 



]T dist 2 (^ e (« e )(0 - . z) < 

(WeOt a -)J. 
4tt 2 4 

~2~7 



]T dist 2 ^9 £ (v £ )(l) - ^9 £ {v £ ){k) , Z) 

< (l + o(l))XF e (<; £ ). 



Since ^fi^ + (s^- , s a - ) C (Ag )°- , we are in a position to introduce the sequences 



for all I G ng s , 



AW 



obtaining 



' £ ' I log <S e I"" 1 
S-C^(0<;j, X! dist 2 (^^(« £ )(0-^^(^)(fc),Z 



Therefore, by (45) we deduce that property i) of Definition [7] is satisfied. In 
order to check that also property ii) of Definition|7]holds, it is enough to prove that 
\\fi e — t £ r) £ \\ — > in X for some t £ — > 1. Arguing as in the proof of Xy e -< QL e it is 
enough to check that |loge|' 1 ||u>(e l2mie ) — w(v e )\\\ — > 0; we skip the details, that 
can be easily checked by the reader. □ 



5.3. Proof of QC e X ST> € . Since XY e < 4:ir 2 SD £ pointwise, we immediately de- 
duce that Xy e < ST> £ . Therefore, the desired order relation is obtained by proving 

GC e i xy s . 

We first observe that by Lemma 2 in [3], there exists a constant C > such that, 
for every v £ AXy e 

% [ W(w(v)) < XY £ (v). 

Let A e y 1 be such that X e il c f2 e (we recall that f2 is star-shaped), and such that 
A e > 1 — - for some constant c. Given a sequence t e \ we have 

(46) (1 + t s )XY e (v) > jT ^ (^|V W ( V )| 2 + ^W(w(v))) dx, 
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for every v G AXy e . Let now [i e be a sequence such that Xy e (fi e ) < C, and let v e 
be such that ^^ h XY e (v £ ) - Xy e (fj, e ) -s- 0. Set then <5 e = e, 



u> e (x) := w(v £ (X £ x)) for every x G f2, 



1 Jw £ 
|loge| h_1 7r 



Then by ( 46 ) we get 

(47) XY £ {v £ )> Y ^-GL c e t '{w £ 



From (47 1, we easily deduce that QC £ C -< Xy e . Indeed, Property i) of Definition 
[7] is a direct consequence of ( 47 ) , while the proof of Property ii) follows as in the 



proof of Xy £ -< QC £ . Finally, choosing t £ — > such that s £ :— Ct £ satisfies the 



assumptions of Proposition 14 we conclude that QC £ ~ G£? te di Xy s ^ SV £ . □ 



Remark 18. In [22] it is proved that for h = 1 the functionals QC £ in (20 1 are 



equi-coercive, and T-converge to the functional 7r|/i|(f2). In [30 [Theorem 3.4]) the 



same r-convergence result is proved for the functionals ST> £ defined in ( 22 1 . In view 
of this result, of Theorem [15] and of Theorem [10] we obtain a new proof for the 
compactness of the jacobians given in [22 , and of the corresponding r-convergence 
result of Ginzburg-Landau functionals in the logarithmic regime. 

Remark 19. For latter use we observe that, if we restrict the Ginzburg-Landau 
functionals GL £ to the fields w £ valued in B\, then the equivalence result stated in 
Theorem [15] still holds true. Actually, given a sequence w £ with finite energy we 
can always project it on B\, without increasing its energy, and without changing 
the limiting behaviour of the corresponding topological singularities. 

This L°°-bound will simplify the proof of compactness properties of the quantity 
j £ associated with w £l we will deal with in the next Section. 

6. New results for the asymptotic of SD e 



As explained in the Introduction (see also Remark 18 1, the T-limit of the func- 
tionals ST> £ is known only for h = 1. On the other hand, the analogous result 
for the Ginzburg-Landau functionals QC £ has been proved by Jerrard and Soner 
in [23J for all values of h (h — 1 and h — 2 being the most relevant cases). In 
this section we use the variational equivalence argument to deduce r-convergence 
results for the screw dislocation model in the | loge| 2 scaling regime. We recall that 
this energy scaling has been already considered in the context of interacting edge 
dislocations in [19], providing in the limit a macroscopic strain gradient model for 
plasticity. We will extend this result to our discrete model of screw dislocations 
without imposing, as in |19j . that the minimal distance between the dislocations is 
of order p » e . 

Before giving the rigorous results, let us explain by heuristic arguments why the 
I loge| 2 energetic regime is somehow critical, and hence gives rise to an interesting 
macroscopic limit. Assume that in the crystal there is a distribution \i of a certain 
number N £ of screw dislocations of unit length. The self energy of the system is, in 
first approximation, proportional to N £ \ loge|. On the other hand, each dislocation 
induces also a far field: the macroscopic strain field f3 has to satisfy the kinematic 
constrain curl /? = //. The elastic energy depends quadratically on (3, and then it 
is proportional to N E . We deduce that if N £ sa | loge\ then the self energy and the 
far field energy corresponding to the macroscopic strain (3 (namely the interaction 
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energy) are of the same order | log e | 2 . Therefore, in the limit as e — *■ 0, the energy 
is given by the sum of these two contributions, a self energy, one-homogeneous 
with respect to the dislocation density, and an interaction energy, quadratic with 
respect to the macroscopic fields /3's satisfying the kinematic relation curl (3 = /a. 
This kind of energy can be settled in the recent strain gradient theories for plasticity 
introduced in [18] . 

We underline that, as it will be clear in our analysis (sec Theorem 23 1 the same 
result holds true for the XY spin system model. 



6.1. The T-convergence result for QL e in the |loge| 2 regime. Here we recall 
the T-convergence result for the functionals QC e in the energetic regime GL e w 
| log e | 2 corresponding to h = 2 given by Jerrard and Soner in [23 . For the sake 
of simplicity we will specialize the results assuming that the order parameters w E 
take values in B\ (see Remark 19 1. 
Given w E H 1 (Q,;Bi) 1 set 

j(w) := w x Vw = (wi(w 2 ) Xl - w 2 (wi) Xl ,wi(w 2 ) X2 - J«2(wi)x 2 )- 

Note that by definition we have J(w) = |curl j(w). Consider the functionals 
GL £ : X x L 2 (fl;R 2 ) -> [0,+oo] defined as follows 

(48) GL e (nJ):= 



1 1 log e | 2 | log e| 



' 7r|loge| 

By [531 Theorem 1.1 and 1.2] we deduce the following T-convergence result 

Theorem 20 (Jerrard and Soner, 2002). The functionals GL e : X x i 2 (fi; R 2 ) 
[0,+oo] defined in (48 1 are equi- coercive: if (/i e ,j £ ) is a sequence with bounded 
energy then, up to subsequences, fi £ — > fi for some fx £ X and j e j weakly in 
L 2 (n;R 2 ), for some j S L 2 (n-,R 2 ). 

Moreover, GL e T-converge (with respect to the same topology) to the functional 
GL : X x L 2 (tt; R 2 ) -> [0, +oo], defined as 

(49) GL(a*,j) ■.= n\^m)+ l -J^\j\ 2 , 

if fi is a measure in H~ 1 (fl) and curl j = 2ir[i, and infinity elsewhere. 

From Theorem [20] we immediately deduce the following 

Theorem 21. The functionals QC e : X — > [0,+oo] defined in (20 1 with h = 2 are 
equi-coercive and T-converge, as e — * 0, to the functional QC : X — * [0, +oo] defined 
by 

gZijj) :=7r|A*|(n) + ^inf \j\ 2 dx, j e L 2 (fl;R 2 ) : curl j = 2tt/x}, 
if fi is a measure in _ff _1 (f2), and infinity elsewhere. 
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6.2. New results for homogenizing dislocations in the | loge| 2 regime. From 
the variational equivalence between QC e and ST> e stated in Theorem 15 (see also 



Remark 19 1, we deduce the equivalent result stated in Theorem 21 for the energy 



functionals corresponding to screw dislocations. 

For the reader convenience, we state the T-convergence result for the dislocation 
energy functionals J- E defined in (50 1 according with (22 1 with h = 2, but without 
the pre-factor 4ir 2 that has been useful to compare the SD model with XY and 
GL models, but which has not physical meaning. 

Theorem 22. The functionals T € : X — > [0,+oo], defined by 

1 \K 



(50) 



inf 1 SD e (u), u e ASV 



|loge| 2 

are equi- coercive and T -converge as e 
defined by 

W:=^W(n) + ±irf{/ n 

if fi is a measure in H^ 1 ^), and infinity elsewhere 



= V 



loge 

to the functional T 



X 



[Or 



-oo 



,/3G L 2 (fl;R 2 ), curl (3 = fA, 



Proof. The proof is a straightforward consequence of Theorem 15 and Theorem 
ED □ 



In order to give the analogue of Theorem 20 for the XY and the screw dis- 
locations model, let us associate to any discrete strain (3^, with u £ AST> e , its 
corresponding piecewise constant strain field (3% := P(Ww(u)) in il £ (where P is 
defined component by component as in ([9])), and extend it to zero in Q\Q £ . More- 
over, given v G AXy e , we set j v := 2nP(X7w(^9(v))) in Q E , so that j v = 2irj3g, v y 
and extend it to zero in Q \ Q s . 

We are in a position to introduce the functionals XY e : X x L 2 (Q; W 2 ) — > [0, +00] 
defined as 



(51) XY £ ( M) j) 



inf ■ 



1 



;XY £ (v), v g Axy e (n) 



Jv 



.7. 



Jloge| 2 "" ev "" ' ' I log e| J ' \loge\ 

and the functionals SD e : X x L 2 (fl; M 2 ) -> [0, +00] defined as 



(52) SD e (/i,j) 



u |loge| 2 I log e I I logs] 



The following Theorem establishes the variational equivalence for the functionals 
GL e , XY e and §D e with respect to the strong topology in X and the weak topology 
in L 2 (n;R 2 ) -> R. 

Theorem 23. The functionals GL e , XY e and SD e are variationally equivalent. 



Proof. The proof follows the lines of Theorem 15 One has only to check that for 
each change of variables involved in the proof, the corresponding fields j(w e ), j Ve 
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and 27r/3^ , rescaled by | loge| share the same weak limit in L 2 . Let us check it only 
for the order relation XY £ < GL e , the other order relations being analogous. 

To this purpose, let (/x e ,j e ) be such that GL 6 (/i £ , j e ) < C, and let w £ be such 
that 



J(w e ) 
7r|loge| 



Me 



I loge| 



1 



I log e| 



;GL(w £ ) — Gh £ (/j, £ ,j £ ) — > as s — > 0. 



Let now <5 e and w e be as in the proof of the order relation Xy £ < QC £ in Theorem 
[l5] (with h = 2), satisfying 



XY Ss (v e )<GL s (w e ) + o(l), \\w(v e ) - w e \\ 2 -» 0, 

In order to conclude, it is left to prove that 

1 



1 



|log<5 £ 



Me 



0. 



(53) 

We first prove that 
(54) 

Since we have 
1 



I logs] 



(jK)-^J^Oin L 2 (fl;R 2 ). 



|logej 



(3'W-iW«e))^0 inL 2 (ft,M 2 ). 



|loge| 



C?O e ) - j(w(v e ))) = (w £ - w(v e )) x 



w(^ £ ) 

loge| ^ |loge| 



x V(w e - iu(w e )), 



by Holder inequality and by integration by parts it is easy to deduce ( 54 1 . 

jTo^lOWO) -3v.) ^ in 



Now, to obtain (53 1 it remains to check that 1 
L 2 (il;M. 2 ). To this purpose, set p £ := 5 £ | log S e \ 2 , and let Qs e ,p E be the union of 
<5 e -squares Qi in Cls E such that the oscillation of w(v £ ) on Qi is bounded by p £ . 
Since XYs F {v £ ) < C | log 5 e | 2 , it easily follows that 



l«\ I 



< 



C 

llogXP 



as e — > 0. 



Therefore, since , lo ^ , j(w(v e )) and u ^ s i j Ve are bounded in L 2 (fi;]R 2 ), we get 



(55) 



I log 



in L 2 (fi;M 2 ) as e -> 0. 



To conclude the proof of ( 53 1 it remains to show that 
(56) lim ||j(w(w e ) - j^lUoo^ R 2) 



= 0. 



To this purpose, notice that, on each i5 E -square Qi € ^<5 £ ,p e we have (for e small 
enough) j Vs = Vu>(#(u £ )). Therefore by Taylor expansion 

\j(w(v e ))-j Ve \ = \j(w(v e ))-Vw(9(v E ))\ <C P f = CJ e |log 5 £ | 4 -> 



on Qs c , Pe , that clearly implies (56 1, and this concludes the proof of (53 1. 

Essentially, the same arguments can be used to prove all the order relations 
between the functionals Gh £ , XY £ and §D e , so that we prefer to skip the details. □ 
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From Theorem [23] we immediately deduce the analogous T-convergence result 
for the screw dislocation functionals. In particular, we obtain a new result in the 
context of homogenizing dislocations, generalizing the results in |19j to the scalar 
case of discrete interacting screw dislocations. For the reader convenience, we 



state the result introducing the dislocation energy functionals Q e , defined in (57) 



according with (52 1, but without the pre-factor 2ir 2 in the strain, and without the 
prefactor An 2 . 

Theorem 24. The functionals Q e : X x L 2 (tt;R 2 ) -> [0, +oo] defined by 
(57) g £ (fi,/3):= 

inf\ T -^—^SD s (u),ueASV £ (Q): T -^-=P, Mu ' 



Jlogep '" "' ' " |loge| | log ej J 

are equi-coercive and T-converge, with respect to the strong convergence in X and 
the weak convergence in L 2 (fl;M. 2 ), to the functional Q : X x L 2 (f2;R 2 ) — > [0, +oo] 
defined by 

(58) g(v):=^m)+ 1 -Jjf3\ 2 , 

if fi is a measure in H (p.) and curl /3 = \i, and +oo elsewhere. 



The T- limit Q in ( 58 ) represents the macroscopic energy corresponding to a den- 
sity of dislocations /i and a macroscopic strain [3. The first term in Q represent 
the self energy of the dislocation density, it is 1-homogeneous, and therefore it is 
compatible with concentration of dislocations. Notice that the compatibility condi- 
tion n = curl j3 agrees with concentration of the density of dislocations, whenever 
/i G _ff _1 (fi). This is the case of concentration on lines, according with the well 
known configuration of wall dislocations. The second term is the elastic energy of 
the macroscopic strain /3, and represents an interaction energy of the dislocations. 
This kind of macroscopic energies, depending also on the derivatives of the strain 
f3, are usually referred to as strain gradient theories in plasticity and they have 
been introduced in |18j . The energy Q derived by T-convergence represents then 
a justification of such phcnomcnological theory, and provides explicit self and the 
interaction energy densities. 



7. Further extensions and conclusions 

In this paper we have investigated the variational equivalence of some model 
characterized by the presence of topological singularities. It is natural to ask if our 
method can be extended to other contexts in the huge field of modelling singularities 
and in particular dislocations. 

In this Section we propose some extension of our approach. We first consider the 
so called core radius approach to dislocations in the two dimensional setting and 
then we investigate the case of three dimensional dislocations. 

7.1. The core radius approach. In order to deal with the singularity of the stress 
field around a dislocation, a very fruitful approach consists in removing a region 
of size e around each dislocation, usually referred to as core region, obtaining in 
this way a L 2 integrable field on the complementary domain. This method has 
been exploited very recently in variational models for dislocations [T3], [3D], US]- 
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Its feature is that the discreteness of the problem is carried by the length-scale e, 
representing the atomic distance, while the mathematical framework is continuous. 

Within the core radius approach, many mathematical details have to be fixed in 
order to make the T-convergence problem well posed and doable: for instance in 
|30j it has been introduced a small penalization for the number of e-disks removed 
by the domain; such a penalization plays the role of the potential term W in the 
Ginzburg-Landau energy, ensuring that the number of dislocations is bounded by 
| log e | , and in particular that the measure of the core region tends to zero as e — > 0. 

This approach, as proposed by Bethuel, Brezis and Hclcin, can be very fruitful 
also as a variant of the Ginzburg-Landau approach to vortices. The variational 
equivalence argument introduced in this paper seems to be a natural tool to compare 
the core radius approach to dislocations and vortices with pure discrete approaches, 
that we believe to be equivalent. 

Finally, we aim to comment that the core radius approach has been proposed 
in [7] in order to compute the renormalized energy, i.e., the lower order term in 
the energy of minimizers of the Ginzburg-Landau functionals in the logarithmic 
regime. This inspired the work in |14J . where the authors compute the Taylor 
expansion of the clastic energy of edge dislocations in a plane. The first term in 
the Taylor expansion is the self energy, that for screw dislocations is given (up to 



a pre-factor) by (see Remark ( IS I ) . The lower order term, corresponding 



to the renormalized energy, depends on the mutual distance of the dislocations, 
and models the Peach-Kohler attractive and repulsive forces between dislocations. 
In our opinion it would be very interesting to investigate the equivalence of this 
renormalized energy for vortices and dislocations within a L-convergence analysis. 
In order to do that, it seems necessary to exploit the equivalence between vortices 
and dislocations for the first order term in the L-limit expansion of the energy 
functionals. Indeed, Definition[7]could be generalized to compare lower order terms 
in the energy, in the spirit of the theory introduced by Braides and Truvskinovsky 
in \1'2\ . The analysis required to compare the renormalized energy for vortices and 
dislocations seems to be very challenging. 

7.2. Three dimensional dislocations. Screw dislocations are essentially straight 
dislocations lines with parallel Burgers vector. The general case accounts more 
complexity, dealing with general closed loops of dislocations with given Burgers 
vectors. 

Here we aim to introduce the elastic energy in this three dimensional setting, and 
the corresponding Ginzburg-Landau energy functionals, formally obtained arguing 
in analogy with what we have done in the anti-planar setting. We will not address 
the problem of proving rigorous equivalence results between elastic energies and 
Ginzburg-Landau energies in this three dimensional context. 

We consider the basic case of a cubic lattice, whose elasticity tensor we denote 
by C £ R 9 x M 9 . Let C K 3 be the reference configuration of the three dimensional 
crystal. The displacement is now a vectorial function from :— il n si? to M 3 , 
and, following the approach in [5], the dislocations are defined on the class of 
two-cells of the crystal. In this three dimensional case, we have three kind of two- 
cells, corresponding to the three different slip planes of the cubic lattice. A basic 
dislocation loop is identified with a pair, given by a two cell and a Burgers vector 
b, that in the cubic case belongs to the canonical base of R 3 . We have then three 
kind of topological singularities, corresponding to the three different slip planes, 
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and each kind of singularity has a vector nature, being associated to a Burgers 
vector in M 3 . It is then natural to set up a Ginzburg-Landau model for dislocations 
considering the vector valued maps taking values near a three dimensional torus. 

More precisely, consider the three dimensional torus T := M. 3 /(eZ 3 ), i.e., M 3 
where we identify x and y if and only if x — y G eZ , and consider the fields 
w : £1 — ► T x [0, +oo). We denote by u £ T the first three components of w, playing 
the role of angular components of a continuous field in classical Ginzburg-Landau 
theories, and representing in our model the displacement function, and we denote 
by p the fourth one, representing the radial component of w. Therefore, for every 
w £ i? 1 (f2,T x [0, +oo)) we define the Ginzburg-Landau energy functionals by 



In our opinion such functionals provide a good material dependent Ginzburg- 
Landau model for dislocations, and this could be justified by showing that they 
are indeed equivalent to suitable discrete elastic energy functionals, defined for in- 
stance according to the formalism introduced in [5] . The rigorous formalization and 
proof of such a statement would require a specific analysis, that is not addressed 
in this paper. Clearly, more general crystal lattices could be considered within this 
approach. 

Exploiting this Ginzburg-Landau approach to dislocations provides a motivation 
to investigate the asymptotic behaviour of these Ginzburg-Landau functionals as 
e — > 0, generalizing in this way the analysis done in [T] , [7] , [22] , [23], [32], [33J to the 
case of vector valued fields whose singularities belong to a given group, and whose 
energy is not isotropic, and not even coercive. 
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